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ABSTRACT 


This thesis theoretically presents the phenomena in- 
volved in the flow of an incomoressibdle fluid over a wedge 
with a second incomoressible, lighter, and less viscous 
Fluid blown through the surface of the wedge. A method is 
developed to determine the inner fluid layer thickness, the 
wall shear stress and the resulting local drag reduction. 


The results predict substantial drag reduction. 
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c COnsStant 

En internal energy 

£59) dimensionless stream function inner-fluid 

£, (nz) dimensionless stream function outer-fluid 

g (x) scale factor for y 

K integration constant 

L reference length 

m constant defined by = 

m mass flow rate 

D local pressure 

oR Prandtl number 

Pomac stagnation pressure 

Q heat into the system _ : 
R fluid properties ratio defined by R = ai 
Re Reynolds number 
Rg gas constant 

D local temperature 

es wall temnperature 

i temperature in the free stream 

u velocity component tangential to surface 

eC xX ) potential speed 


U reference speed 
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SHY 
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VelocLty Cemponent normal to surface 
blowing velocity at the wall 
velocity component normal to u and v 
velocity 

work done by the system 

arc length along the surface 

body force in x-direction 

direction normal to surface 

dummy variable for f 

dummy variable for ¢' 

dummy variable for f'' 

body force in y-direction 

direction normal to x and y 

body force in z-direction 


L : d 
SS & Glxg (xy 


Similarity constant given by 
"dimensionless" wedge angle (angle is 87} 
boundary layer thickness 

dimensionless similarity variable 
dimensionless inner fluid laver thickness 
dimensionless parameter 

mass density 

dynamic viscosity 

kinematic viscosity 


stream function 


shear stress 








shear stress at the wall 


: reference shear stress 
oo 








ACKNOWLEDGMENT 


The author expresses his appreciation to Dr. J. V. 
Sanders for leading him through the darkness of scientific 
research and for his availability to assist in breaking 
meaglocks, and to Dr. F. J. Jouaillec who, by his fresh 
ideas in these trying times, was able to give impulses to 
the progress. 

Special thanks are given to my wife Carla van Dord- 
zeegers and our three sons Bas, Lex and Rob for their 


patience and endurance. 





hee oe ROD UCT ON 


The ever increasing operating depths and speeds of modern 
submarines and the much higher resistance to the effects of 
contact or nearby underwater explosions all require a dra- 
matic increase in the speed of underwater weapons. The 
existing power plant systems seem to have reached their 
limits in power density and radical changes in body shapes 
are limited by launching facilities. The answer might lie 
in the field of boundary layer control. Possibilities are 
injection of polymers into the boundary layer to maintain 
laminar flow, suction to remove votential turbulence, more 
sophisticated weapon shapes to reduce form drag, or intro- 
ducing a gas layer between the surface of the body and the 
irauid. 

Of all the options, the two-phase boundary layer seems 
to be the most promising possibility for obtaining an order 
of magnitude increase in speed. 

In this thesis, a method is developed to calculate the 
properties of the laminar, two-phase flow produced by in- 
meeting a light, incompressible low-viscosity fluid into the 
boundary layer of the incompressible flow over a wedge. This 
is just a step on the long way to a high speed underwater 


weapon. 


Ne, 





If. BACKGROUND 


A. DISCUSSION OF RELATED PAPERS 

An extensive literature search was carried out and the 
relevant papers are discussed below. The reduction of drag 
by use of polymers is deleted from this survey because the 
drag reduction expected by this method is probably not suf- 
ficient to increase the maximum speed of the underwater 
vehicle substantially. Amazingly enough, the most recent 
paper in the field of two-fluid laminar boundary layer flow 
was published back in 1969. 

mod. Pretsch [1] 

Pretsch starts with Prandtl's boundary layer equa- 


tions for single-phase flow with suction or blowing: 


du du aU o°u 
ax 








and boundary conditions y = 0:u2=0, v = Vee yvRe>~> :; u=0 
where u is the velocity tangential to the surface, v is the 
merocity normal to the surface, x is the arc length along 
Mie surface, y is the normal to the surface, U(x) is the po- 
tential velocity, v is kinematic viscosity, Verte) is the 
normal velocity or blowing velocity at the surface and Re is 


the Reynolds number, Re = a 


al 





Hommceronagusictaon (vy. + —~) the above equations 


0 


Simplify to a second order, ordinary differential equation 


of the form 





Z yv 
ole du _ See W 
a ate ae = 0 where e 
ae Vv 
and boundary conditions §— = 0: u = 0 


ee = (x) 


The solution of this equation is the velocity profile 





For strong blowing we > ©) the assumption is made 
that the viscosity term can be neglected and the above equa- 
tions again simplify to a second order, ordinary differential 
equation but now nonlinear. 

For a number of values for similarity parameters, 
Pretsch gives analytical expressions for solutions of flat 
plate and wedge flow. 

The assumption of negligible viscosity is too re- 
strictive for the case under investigation in this thesis. 

wae. oJ. Natson [2] 

This extensive report consists of three parts, i.e. 
Theory of Similar Velocity Distributions, Flow with Uniform 
Merton, and Flow with Variation of Suction Velocity. The 


report is based on the asymptotic theory of boundary layer 
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flow. Aithough the author is aware of the existence of 
Pretsch's work, no comparison of results is made. 

In Part I the author arrives at a number of tables 
giving velocity distribution and coefficients of the series 
expansions as functions of similarity related constants and 
dimensionless parameters. 

In Part II the similarity principle is abandoned but 
the suction velocity is uniform. The velocity 1s expanded 
into a series from which series for momentum and displace- 
ment thickness and skin friction are derived. In the sub- 
sequent section a number of applications are investigated 
and the results are totalized in a number of tables. 

In Part III the asymptotic theory is extended to 
cover general two dimensional flow with arbitrary distribu- 
tions of mainstream velocity and suction velocity. The 
results are again displayed in tables. The author also de- 
rives the results in Part I as a special case of Part III. 

As this report only deals with suction it is not of 
direct interest to this study. 

Peek. OD. Cess and E. M. Sparrow [3] 

This paper deals with a two-phase flow in forced 
convection over a flat plate. The inner fluid is obtained 
by boiling the outer fluid at the wall so that a vapor film 
is formed. The authors arrive at the following conservation 


equations: 
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Bee eee! CCS. 0 
fg Gg 
STE 
eG where 9 = & Fr 
W ee) 
tS + foth = 0 


with boundary conditions: 


il 
© 
CD 
>) 

ll 
t— 


plate surface: eres (0) 


g 
5 
(ou) 
i) ee) ah) 
U G Gl 


interface: £ (0) = 7 - 
1 =— —- 1 
tS, (0) “g (Ngi) 
lL 
(ou) ) 
tt = an | tt 
f. (0) J fy (Ni) 
9 (n 3} = 0 
free stream: to -~ 2as Ny, se 2 


where EG and fF are the stream functions for respectively 
inner and outer fluid, T 1s the local temperature, ie is the 
Semeperature at the wall, T._ 1s the temperature in the free 


stream, PR is the Prandtl number, 0 is the mass density and 


fers tne dynamic viscosity. ra is a function of the dimen- 
rU 
= Be Sioa : 
Sionless parameter Ng defined by Ng 24] Tg : Nei is the 


value for this parameter at the interface of the two fluids. 


f. is a function of the dimensionless parameter n, defined 


L 

U 

by 1 = a The subscripts L and g refer to respective- 
L 

ly outer and inner fluid. The primes denote differentiation 


with respect to the relevant dimensionless parameter. This 
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oroblem is solved by using an Emmons-Leigh [4] table anda 
successive approximation procedure. 
4. K. Nickel [5] 

iimewissoape., the laminar flow of an incompressible 
fluid without heat transfer is investigated using a method 
of bounding the solutions applied to Prandtl's equations. 
After introducing the well-known dimensionless parameters 17 
and f(n), Prandtl's partial differential equations reduce to 


the ordinary, third order non-linear differential equation: 


ieee varies fio) =o 


with boundary conditions: £(0) = Eo 


tee (LO) 
ae 2) 

The solutions are reflected in a number of graphs, 
especially "“abbildung 8" giving £''(0) as a function of 8 
and a will prove to be very useful in this thesis. 

See Rk. D. Cess and E. M. Sparrow [6] 





This paper is an extension of the paper discussed in 
Section 3 above and deals with the case of a subcooled fluid. 
To extend the earlier results 1t was necessary to consider 
the conservation of energy within the fluid. The results are 
again reflected in a number of figures. 
ese wie oDaGnOw, V. K. Jonsson and E. R. G. Eckert [7] 
In this paper the laminar two-phase flow, obtained by 


injecting gas into a liquid at the wall of a flat plate, was 


IES: 





Mavyestigqazted. The authors arrive at the following ordinary 


differential equations 


tt 
oO 


: os £ ft 
mor the gas g gtg 


tt 
Oo 


mere the Liquid: fF + fot. 


with boundary conditions: 


at the wall: £, (0) =f (blowing parameter) 


gO 
25, = 0 (no slip condition) 


at the interface: Ey ing3) = £, (90) = 0 (no mass transfer) 
EG (N43) = R fT. (0) (stress continuous) 
1. 
(ow), : 
where R = \--~— 
(pu) 
g 
t. (154 = f. (0) (velocity continuous) 
Obg i (C FL (ny) pECuVGlOcLUy approaches U_) 


The authors solve these equations by applying a 
series expansion and subsequently a successive approximation 
procedure, similar to the one described in the paper by Cess 
and Sparrow (Section 3 above). 

The results are given ina table and a graph. It 1s 
proven that evaporation at the interface only works 
advantageous. 

fee. Pachenberger and J. D. Offutt [3] 

Under assumption of 

- shear stress in liguid is negligible 


- interfacial velocity and pressure distribution are those 
imposed by the potential flow 
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- the inertia of the gas is negligible 
- the gas flow is laminar 
- body moves steadily through water at zero incidence 


- gas and water properties are constant; in particular the 
density of the gas remains constant 


the problem of a two-phase flow over different kinds of 
bodies is solved by a finite differences method. Results 
are given in a number of graphs. 
Pei mor Scere Ono. Sabkdollsand J. T. Byrne [9] 
This paper discusses experiments with a "gassing" 
model to investigate possible drag reductions. The vapor 
layers were obtained by film boiling, sublimation of a solid 
surface and chemical reactions. The exveriments indicated 
that for film boiling a region exists for which a stable 
vapor layer may be expected. 
Assuming 
- steady, two dimensional laminar flow 
- pressure everywhere constant 
~ heated surface temperature constant 


- temperature of liquid vapor interface is equal to satu- 
rated liquid temperature 


- thermodynamic and transport properties of liquid and 
Vapor are constant 


eeemermal radiation effect on yarOr and liquid layers is 
neglected 


- buoyancy effect on forced convection boundary layers is 
neglected 


-~ velocity and temperature profiles are linear in vapor and 
parabolic Thigts Ar alke Rulalys 


1 





and using the Pohlhausen integral method the results of the 
experiments are predicted. In some of the experiments dra- 
matic instabilities were found, and the authors had to con- 
clude that producing a stable vapor layer between a surface 
and a liquid drastically reduces friction drag but unstable, 
two-phase flow (like nucleate boiling and uncontrolled sub- 
limation) may actually cause the drag to increase. 
cele. ang [10] 

In this patent Lang "designed" a torpedo capable of 
using the gas film drag reduction phenomena. The lift is 
provided by the body shape and the exhaust gases are used to 
produce the gas layer. At the rear end of the torpedo the 
gas layer is removed by suction. 

Moe &. R. van Driest [11] 

The author indicates that the quest for higher 
speeds of underwater craft can be solved by using the knowl- 
edge acquired in high speed aerodynamics. 

In a number of figures the author gives local fric- 
tion coefficients as a function of Reynolds number, speed, 
wall temperature and ocean depth. Further effects of hydro- 
dynamic heating, cavitation and insertion of additives are 
discussed. 

The author concludes that although the propulsion 
systems have to be developed, it is, in principle, possible 


to obtain high velocities for underwater vehicles. 
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ee QUO VADIS 

Considering the above discussion it was felt that more 
knowledge had to be obtained of the phenomena taking place 
in the two-phase boundary layer. A first step would be to 
extend the work on flat plates to two-dimensional bodies 
with pressure-gradients. It is logical to follow the paper 
of Sparrow, Jonsson and Eckert [7], discussed above in Sec- 


tion A-6. 


LS 





eee, 


A. GOVERNING EQUATIONS 


THEORY 


To predict the skin friction and drag characterics from 


first principles the starting position is based on the basic 


conservation laws governing the velocity distribution, here 


in Cartesian coordinates, 


ae 30 3 {0u) 
Sentinulty: rs te 25% te 
Bue Boe 
Navier-Stokes: rs = X ae ~ 








a 4 
oC 
zs 

Ben of State: D= oO Ry AB 
= ate, T) 


Bpercy Equation: =: = = + = 














elie 
NIG 
a 


xe . e s e 
where w 1s the velocity with components u, v, w in respec- 


tively the x, y and z direction, 
e) 3 J 0 
fae ox * “ay * “Sz 


20 


po 
Dt 


is the operator 


, &, Y and Z are the body forces in 





Pee cemvelyecwey ana 2 directions, P is the local pressure, 


is 


R. is the gas constant, Q 1s the heat into the system, Ee 


the internal energy and W is the work done by the system. 
Assuming steady state, absence of body forces, viscosity 

oP . 2 ee 

ay ms negligible, 


laminar flow, and both fluids incompressible (the energy 


u(T,P) is constant, two-dimensional flow, 


equation decouples), these equations simplify to the Prandtl 


boundary layer equations: 


du OV 
@Sencinuity: aot ay * 0 
du du 1 93P a 
Navier-Stokes: u-—— + v— =- == + ws, 
ox dy a) 8 ors 34° 


Using the subscript g for properties of the inner fluid and 
me@emsunoscript L for the outer fluid, it is possible to write 
two sets of equations describing the two-fluid laminar 


boundary layer. 








Inner fluid: —— + —= 0 (1) 
ox JV 
9U du 1 3P 56 
an + V=—=- aes < V oy (2) 
oy e) Jay 
Serer fluid: ey ee 0 (3) 
ax 3V 
you ~ you = = a oP + VY ar (4) 
opie JY oo x Ley 


To describe the problem fully, however, it is required to 
state the boundary conditions and the conditions at the 


interface between the two fluids. At the wall, in the inner 


a 





mimid, the velocity parallel to the wall, u, 1s zero (no 

Blin) and the velocity perpendicular to the wali is equal to 
me imposed blowing velocity. At the interface the tangen- 
tial velocity and shear stress have to be continuous. Fur- 
thermore, the mass transfer across the interface is assumed 
mene Zero. In the outer fluid u has to approach U(x) as y 


approaches infinity. A formal statement of these conditions 














is: 
at y = 0 u = 0 

V= Vv 

W 

at interface uy = uy 

' be = JU, 

gi oy Li oy); 
m= 0or Vg (interface) = Vy (interface) = 0 


Memourer fluid u- U(x) as y > » 

For the given assumptions, equations 1-4 and the boundary 
conditions describe the problem completely. This system pre- 
sents substantial difficulties in attempting a general solu- 


tion. Similarity solutions, however, will solve the deadlock. 


See LRANSFORMED EQUATIONS 

Similarity-type solutions are well known and have the 
characteristic that the velocity profiles have similar shapes 
at all x-positions along the surface. Mathematically this 


gives the enormous simplification that the governing partial 


ee 





differential equations are converted into ordinary differen- 
tial equations. Before defining the dimensionless para- 
meters it should be noted that the continuity equations are 


satisfied by the stream-function y» defined by 


GS 
ll 
ales 
~ FS 
~ 
3 
Qu 
< 
iH} 
{ 


Following Schlichting [12], the following dimensionless 


parameters can be introduced for the inner fluid: 


UL 
Reynolds number RE, iiervane (where U, and L are some 
g 
reference quantities) 
y/Re 
Ney = aes Eg = -- (where g(x) is a dimensionless scale 


PacwOm LOr y) 


w(x,y) /Re 
E(n,,é,) = —————_4 or wlx,y) = £(n 


c : yi U(x) g(x) 
g ee One) (35) % 


g aT 
Re 
eG 


t 


g 


memsimplitcy further manipulation, it is wise to state at 
this point that to ensure similarity = ee can not depend 
Smee = & (x) so that £ = Ein). 


g g 
From the above, the velocity components may be derived, 


u = U(x) £ ‘a 
Tee i) 
and v=- = =~ (U(x) g(x)) + GU es es) n ac Ng ) Ee = g(x) 


g 
—— / Re, 
Inserting these expressions into the Prandtl boundary layer 


equation and after a considerable amount of algebra the equa- 


mon Converts into 





2. 
1d 1? - Lia f t 
fy +af ft, B ( g ) 


li 
© 


where a = 2g) (U(x) g(x) ) 


g = L g? (x) au (x) 
U 


Gise 


Reasoning in the same way for the outer fluid after 
defining 


= ( 


= Mae) fe - = x 
L Vv a oes) (Poe ae cl L 
i - p(x,y) / Rey 
i'"n75,) = Dx) g(x) 
one arrives at 


ee 14 = 12 
fT. + af, ft. ee fe. ) 


i 
= 


where a and 8 are defined the same as for the inner fluid. 
It should be noted that to satisfy the similarity re- 

quirement : and 8 should not depend on x. Still following 

Seer chting [12], it is possible to arrive at the following 


@emeral conditions for U(x) and g(x): 








8 
LN 
U (x) Dae | x 2a-6 ; 
5 = K | (2c - 3)F where K 1s an integration constant 
x 2 
gix) = _[ (2 SS) > Ux) 
Now, without loss of generality, take «a = 1 and introduce a 
new constant m = —- Soee lot eeneaexpresslOns LOr cts) and 


oO 


ax) convert to 


24 











Ties), Glee x)" 
She ir) 1) 
(ee ee 
J ~ Afm+l L U(x) 


Furthermore, the transformation equations for the ordinates 


= mt+l U(x) 2 m+1 Z a 
are Ng = y PPE and ny = (y-d) Go 





VX 
Z 


With Schlichting [12], it can be noted that these re- 
quirements can be satisfied by the flow in the neighbourhood 
of the stagnation point of a wedge with included angle 37, 
where U(x) 1S given by cx 

Next the boundary and interface conditions have to be 
converted by expressing them into the newly defined dimen- 


Sionless parameters 





At y = 0 
= 0 = U ics = 0 : = 
or Ng 1b (x) g (n) or ty (QO) 0 
= ee Sige =A): 
Ve My = ZF SF u(x) g(x)) 
¥ Re 
-v_.. /Re 
f EVE ware 
oF Fg(ng) = “ate gad 
abt 
at interface g = Nei MoH 70 eor Ey (gi) = FL (0) = 0 
nN, = 0 
Uy = Uy or 25 (N34) = FL (0) 
du du, (ou), bs 
=o a! f fi )s Foot : = 
Hg dy Uy, ay Or g (n,) R re (0) where R (OH) 


ig ~ or) Teeny yn +o 
fey - Or n+ deo Ut) Or a 7 aS Tr, , 


aS 





It should be noted that 


- to ensure that £5 (Mg) is independent of x: 


nN He (x) ) : Inserting the expressions for U(x) 





and g(x) one obtains vee, OTs) 


- R depends on the properties of both fluids depending on 
temperature and pressure. 


Summarizing the problem can be formally stated by 


2 
fc ee ea 1- £_' = 0 5 
equations 5 5-6 B ( g ) (5a) 
2 
££ ttt tr = ' = 
fF - fs. fs. + &(] EG ) 0 (5b) 
mie COnaditions: (0) = 0 
a _ U (x) 
£{0) Yo) URL) Vv ane where vv Fania 
g 
— — 
eG "gi! (9) =~ 0 
Bo Uitga) SH O) 
Es (N64) = R tS (8), 
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Vee exe bDOreED WAYS OF SOLUTION 


Pee OERIES SOLUTION 
Following the procedure of Sparrow, Jonsson and Eckert 
[7], a series solution of the momentum equation seems to be 


appropriate. Expanding ale, with respect to ‘ls = 0 gives 


See (Cem: aun O) ne 
at: elo 


= 8 + + + 
£ (N14) E (0) + Se ( ie aT an 


The first two terms in the expansion are available from the 
boundary conditions. The eo) is unknown for now and will 
be determined by the solution. Higher derivatives can be 
expressed in terms of lower derivatives by taking sequen- 


tially derivatives of both sides of the momentum equation 


€''! = -£ ft! = 3(1 - £'7) 
and evaluating at Ng = Q. 
Some expressions for derivatives are: 
aoe oe= (oy (oy = 

e'FV) (9) = £7 (0)£'"(0) + 3£(0) 

Moe (VEL (0) = @£-(0) + (26-1) £'*7 (0) 
eT) (9) = £7(0)£''(0) + 8£°(0) - (86-5)£(0)£''7(0) 

- (68-4) 8£'" (0) 


re | 





meawas impossible to find a recurrency in the expressions. 


On applying the interface condition £(n 3) = 0 makes it pos- 


sible to write 


=f (0) 
g a VS) Tig) e ge 


1 fo CO ne! 
g : g (QO) 


nw 
li 


It is now possible to express 2g 1 in liguid and gas para- 
meters. Here Sparrow, Jonsson and Eckert refer to an earlier 
paper by Cess and Sparrow [3]. This paper is discussed in 
Chapter II, Section A3. This paper gives a table relating 
ff,’ (0) and f° (0) which corresponds to a solution of the 
Semactlon for a flat plat. 

At this point it was realized that the data necessary to 
go any further was not available, therefore this method was 


abandoned, and other methods of solution were explored. 


B. NUMERICAL SOLUTIONS 
mee LODE [13] 
ue omen eeiactLVe Ordinary differential equation 
solver developved by the Computer Center of the Naval Post- 
graduate School. The program accepts a number of first order 
non-linear ordinary differential equations with associated 
initial conditions and solves the equations. 


In this case equations 5 were converted into: 
——— = ' ° = 
ae Y>5 ie Note: y, f (6a) 
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dy. 





a a i (6b) 
Bike 2 
a ee GR eee, ot 


where ie ee and y3 are dummy variables for the derivatives 


Memeetn), With initial conditions 


y,(0) = £(0) (7a) 
y,(0) = £'(0) (7b) 
y,(0) = £''(0) (7c) 


The equations 6 are the same for both fluids but as noted 
before the initial conditions differ. 

Sanders [14] suggested that it should be possible to 
guess the thickness of the inner fluid layer Nei and then a 


value of ee could be found by trial and error so that 


no mass transfer occurs across the interface (EM) = Q). 
The values for te —c and i's at the interface give the 
gnitial conditions for the outer fluid equations 6. The 


values of zr, (My) should approach unity for large enough Mn 
while Enz) approaches 0. If tM yz) does not go to unity 
the initial guess for the layer thickness was wrong and 
another value for the inner fluid layer thickness Teh 1s 
chosen. This process is repeated until a layer thickness is 
Beund Lor which ff) = 1 and fff) = 0. 


mae Calculations Carried out on the IBM 3033AP com- 


puter system showed that this procedure worked, but while the 


PLS, 






calculations were fast the interactive side of the program 
has proven to be too time consuming. 
2. Desk Top Computer 

So the necessity arose to find some way to decrease 
the time spent communicating with the computer and to auto- 
mate the trial and error process. Furthermore, it was 
necessary to find some way to interruot the execution and 
to execute the program by steps to determine instabilities 
in the solution. 

Since a Wang 2200 desk-top computer with line printer 
was available, the problem was rewritten for this computer 
and the solution for a number of cases was carried out. It 
has to be realized that a desk-ton computer is relatively 
slow so that a compromise between calculation time and ac- 
curacy of the result had to be found. It was decided that 
the primary goal was to find a method of solution and not to 
moaiate a large amount of accurate data. In the light of 


all the assumptions high accuracy would be doubtful anyway. 
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V. SOLUTIONS 


Pee DEFINITION OF THE PROBLEM 
The momentum equation will be written as three first 


order differential equations as for the "IODE program": 


dy, 
Tao. 
dy 

Dike = Pei 
mes oS 
dy 

ae et 2) 
een ee wa 


Now multiplying both sides of these equations by dn and 
changing to finite differences one arrives at 
Ay, = y5°An 
AY 5 = y3°4n 
Vee Nae as wee 8 *) an 
3 1°¥3 ea 
Now by stepping n through the desired range of values, 


g 
and Y3 it is possible to calculate Yy" Yo and Y3 as a fone — 


Seeeting with known or guessed initial values for n_, Noe aS 


mon Of n. 

It has to be noted at this point that we approximate the 
integrals by a "Lower Riemann Sum". It is known, however, 
that using a "sufficient" number of steps guarantees a satis- 
Bactory accuracy. This procedure is applicable for both the 


mum@emercanad the outer fluid. The initial values for the outer 


Sh 





med follow logically from the solution of the inner fluid 
and the interface conditions. 

Given values for the blowing parameter £ (0) and the 
Meeperties of the fluids parameter R the procedure is as 


follows: 
- Choose an inner fluid layer thickness. 


- Find £''(0) by an interval halving method so that f£(n 3) 
(no mass transfer across the interface) equals 0, or <I 
in our program is smaller than a chosen value. 


Note: £(0) = given; £'(0) = 0 and f£''(0) is chosen as 
indicated. 


- Apply the interface conditions on the values for f, £' 
and f'' at the interface to find the initial values for 
the solution of the outer fluid equations. Choose an 
end value for n, large enough to be well "outside" the 
boundary layer. 


- Check if Se ee =— Met emensoutter. '(°) = 0 as it is 
possible that an unstable solution is found. In case 
the conditions are satisfied the solution is found. If 


not another guess on the inner fluid layer thickness has 
to be made. 


- The value of the inner fluid layer thickness is chosen 
in another interval halving process. 


- This process is repeated until the requirements of 
£,' (=) = 1 and £,'' (2) = 0 are met. 


A copy of this program is enclosed as Appendix A. 
From the above, it follows that the initial intervals 


© 


KX 


2 and Ned have to be determined. The interval for 
"gi Should be as narrow as possible and the interval for 
£''(0) should be a little wider than expected to prevent in- 
Stabilities and "limit-riding". These values were determined 


by trial and error and the results are discussed ina later 


chapter. 
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eee TESTING THE PROGRAM 
ime Number On Steps Of Independent Variable 
The program is used with 100 steps of the independent 
variable. A sample problem was solved with both 100 and 200 
steps. The deviations were in the fourth decimal (see Table 
5.1). The time required for the calculation however was 
doubled from approximately 45 minutes to 90 minutes. 
2. One Fluid Flow 
It is clear that when R approaches unity, the two 
fluids are hydrodynamically identical. Computer solutions 
were obtained for a single fluid flow for 8 = 0 (flat plate) 
and 8 = 1 (stagnation flow) and further for 8 = 0.2, 0.4, 
0.6 and 0.8. The results were in good agreement with Nickel 
fie abbildung 8". (See Table 5.2.) 
Pepe colution for R Close to l 
A solution with a R of 1.1 was compared to that for 
a one-fluid problem and the results were close enough to 
give further confidence in the correctness of the program to 
G@eamake non-unity values of R. (See Table 5.3.) 
meecoOlution for Small Wedge Angles 
mae solutzon for a 6 of 0.01 was close enough to that 
of the flat plate to give confidence about the capabilities 
of the program to correctly handle non-zero values of 8. 


(See Table 5.4.) 
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fee DEFINITION OF A SAMPLE PROBLEM 

To demonstrate the capabilities of the program, a special 
case was chosen for thorough investigation. For this case, 
"air at standard conditions" and water were chosen as the 
two fluids so that R = 200. For a range of blowing para- 
meters C0); and wedge angles 8, the velocity profile was 
calculated and the layer thickness Noi of the inner fluid 
and the shear stress parameter ie 0) at the wall were 


determined. 


Sela LCKNESS OF THE INNER FLUID LAYER 
Figure 6.1 gives the dimensionless thickness of the inner 


fluid layer n as a function of the blowing parameter £., (0) 


Ga 
for various wedge angles 8. Figure 6.2 is an enlargement of 
these curves for small blowing-parameters. These graphs lead 


to the following observations: 


- For small values of the blowing parameter the curves are 
mesa gnt lines, no, decreases very little with increasing 
wedge angle and i§ directly proportional to £_(0). 


- For larger blowing parameters £,(0), 7 4 increases more 
rapidly with increasing £,(0). “For latge blowing para- 
meters and small angles the thickness of the inner fluid 
layer increases drastically with increased blowing 
parameter. 


- For large values of the blowing parameter, the thickness 


of the inner fluid layer decreases rapidly with increas- 
ing wedge angle 8. 
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@. WALL SHEAR STRESS 

Figures 6.3 and 6.4 give the dimensionless wall shear 
stress Eg Ve) as a function of the blowing parameter £(0) 
and the wedge angle 8. 


- For all values of wedge angle 8, blowing decreases the 
wall shear stress fgi' (0). This reduction is spectacular 
for small values of the blowing parameter; e.g. for both 
the flat plate (8 = 0) and stagnation flow (8 = 1) 
£g' | (0) 1s reduced to half the non-blowing value when 
fg (0) =~ 0.0025. 


~ For all values of the blowing parameter the wall shear 
stress £''(0) increases with increasing wedge angle 8; 
the smaller thickness of the inner fluid layer caused by 
meme higher potential velocity. 


~ For large values of the blowing parameter, the wall shear 
Stress is very small and nearly equal for all wedge 
angles. The values obtained indicate that the wall shear 
stress for large blowing parameters approaches the wall 
shear stress for a wedge in an infinite volume of inner 
fluid. (See Table 6.1.) 

PeeeeveLOCITY PROFILES 

megure 6.5 gives velocity profiles (f£'(n) vs. n) for 


Various blowing parameters for a wedge angle 8 = 0.2 (36°). 


The Nr and 1s nave been expressed in the same relative scale 


y/Re yao} Re | 
by noting that le on (x) and 17, = Smeaica resulting in 
ee. 1 
Any, 5 


The following observations can be made: 


- For the smaller blowing parameters the velocity profiles 
in the inner fluid are nearly straight lines. 


- The stronger the blowing the thicker the layer of the 


gnner fluid causing a larger proportion of the potential 
ferocity to be obtained in the inner fluid layer. For 
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small blowing parameters this proportion of the potential 
velocity increases rapidly with increasing blowing para- 
meter. Figure 6.6 gives the proportion of the potential 
mewocity (Or £ '(n7_.)) aS a function of the blowing para- 
meter eg ec af 


- Particularly interesting is the fact that the velocity 
profiles in the inner fluid and the outer fluid have op- 
posite curvatures. This is comparable to the boundary 
layers associated with the gas-filled cavity attached to 
a body after water entry. The cavity around such a body 
has a stable interface; breakdown occurs only when too 
much gas is left at the rear separation. Although no 
rigorous proof can be given it is felt that this charac- 
teristic of the velocity profile provides stability of 
the interface in two fluid boundary layer flow. 


E. DRAG REDUCTION 

Horcnaracterize the amount of drag reduction, the local 
drag parameter sali obtained for the two fluid flow was 
compared to the local drag parameter obtained without 
blowing. 

To compare these local drags, a minor converson has to 


be made. The shear stress at the wall is 


a {du 


ea leases 
(dy y=0 


In terms of the transformed variables, the expression for T 


becomes 


Tug ey /Feg (ENO) 2 


es (ayy eG 


which can be expressed in the outer fluid parameters 


dg 
(ou) 
_ U (x) =— oe 7 
eeu oo mee 1 (on), et (1) 
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For the outer fluid, without blowing, the reference shear 


stress To Sic 


=, —tbllesles 7 
To =U Egtay / Rey, ee 
f (0)=0 
Now taking the ratio: 
= A: ee 
Ce oo. (0) =e0 
g 
The values of this reference shear stress parameter fi (0) 
are calculated and reflected in Figure 6.7. Since Risa 


known constant, for each value of the shear stress parameter 
mf) it is possible to determine the ratio tT/To- Table 
6.2 gives the values of T/T, as a function of wedge angle 3 
and blowing parameter eg! fOran = 200: 

Inspection of this table reveals that a drag reduction 
is observed for all blowing parameters and wedge angles. It 
is remarkable that for small values of the blowing parameter 
mee ratio T/T) decreases tor increasing wedge angle 3 while 
for large values of the blowing parameter £(0) the ratio 


t/t, increases with increasing wedge angle 8. For all wedge 


0 
angles 8, however, the ratio T/T) decreases with increasing 


blowing strength. 


feeeereJULRED POWER 
The power required to maintain blowing is 


Power = P *°* Vv 
W 


where P is the local pressure given by 


Sy 





3P | 
= +) ——— 
cence | ax | 
where PemaG is the stagnation pressure at the leading edge 
oP _ oa Z 2m-1 
of the wedge and ee °g ace ey = m X SO eel che 
inner fluid is at stagnation pressure P a pressure dif- 


STAG 
dP 


ference of ae at the wall exists so that the fluid flows out 
by itself. The flow rate then can be controlled by the 
blowing hole's size. Of course, the ejected fluid has to be 


replaced requiring it to be brought up to P In practice 


STAG © 
Hees 1S of less interest as a considerable amount of exhaust 
gases of the propulsion systems at high pressures should be 


available. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


It has been shown theoretically that a dramatic drag 
reduction is possible for the infinite wedge by blowing with 
a lighter and less viscous fluid, and that the parameters 
involved can be calculated with relatively simple means. 
Furthermore the developed method gives another way of solving 
the flat plate case and the stagnation flow case. 

This work directly invites a number of follow-on investi- 
gations of phenomena associated with: 

- finite two-dimensional bodies 

- three-dimensional bodies 

- a compressible inner fluid 

- experimental verification of the theory 
- stability of the interface 


all of equal importance. 
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APPENDIX A 


TO SOLVE TWO FLUID LAMINAR FLOW CHARACTERISTICS 


REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
OIM 
OIM 


THIS PROGRAM IS DESIGNED TO CALCULATE F,F’,F" 

AND NGI FOR A TWO-LIQUID-LAMINAR-BOUNDARY-LAYER 
GIVEN A BLOWING-PARAMETER F(0O),A WEDGEANGLE B AND 
A LIGUID-PROPERTIES-CONSTANT R. 


TO USE THE PROGRAM ONE HAS TO GUESS ON AN INITIAL 
INTERVAL FOR THE INNER-FLUID-LAYER-THICKNESS AND 
FOR THE LOCAL SHEARSTRESS-PARAMETER F" (0). THESE 
VALUES CAN BE OBTAINED FROM GRAPHS IN THE THESIS. 
ONE HAS TO TAKE THE INTERVAL FOR THE THICKNESS AS 
SMALL AS POSSIBLE AND THE INTERVAL FOR THE SHEAR- 
STRESS A LITTLE BIT LARGER THAN EXPECTED TO PREVENT 
INSTABILITIES AND LIMITRIDING. 


IN SHORT INGERT VALUES FOR 

= F(O) AS Gili), LOWER LIMIT OF F"(O) AS G& AND 
= eee ctnm | Gt FCO) AS GY IN LIRE 120 

= LOWER LIMIT OF NGI AS H6 AND UPPER LIMIT OF 

- NGI AS H/ IN LINE 130 

= WEDGE ANGLE B AND FLUID PROPERTIES CONSTANT R 
= IN LINE 140 


THE OUTPUT WILL CONSIST OF TWO TABLES GIVING 
FG,FG’,FG",FL,FL’ AND FL® AS A FUNCTION OF THE 
DIMENSIONLES PARAMETERS NG AND NL,PLUS A HEADING 
OF THE MOST IMPORTANT PARAMETERS AND CONSTANTS. 


Serguei. Gl(lOr. Ger toOl), G3(101) 
PO Gnectaeeit tory, tettOl), L3i1ol) 


GO(1)=0: Git1)=-.5: Ge(i1)=0: GE=0: G7=2 
H6=.'95: H7=1.05 
G4=G6: GS=G7 


A=1: 
Z=0: 


B=1:2 R=200 
Z1=0 


H8=HE: H'D=H7 

H= (HE +H) /2 

G8=G6: G3=G7 
G3(1)=(GK+G'3) /2 

IF H=Q THEN 400 
meee * 35 THEN 180 
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PRINT “AFTER 35 ITERATIONS IN GAS UNABLE TO MEET REQ’S" 
Bom .G1¢101),G2(101),G3¢1) 
GOTO 495 
HO=H/ 100 
FOR I=1 TO 100 
Hi=G2(1I)+#HO 
H2=G3(1)*HO 
H3=HO* (B#GS(1)#62(1)-B-A¥G1(1)*#G63(T)) 
Gi(I+1)=G1¢1)+Hi 
G2(I+1)=G2(1) +He 
G3¢(I+1)=G3(1)+H3 
GO(I+1)=GO(I) +HO 
NEXT I 
Z=Z+1 
IF ABS(G1(101)) < .O00000001 THEN 255 
ie Gitio1) > © THEN 250 
GB=G3(1): GOTO 150 
G9=G3(1)!: GOTO 150 
LO(1)=O: L1i(41)=Gi(101):L2¢1)=G2(101): L3B(1)=G3(101)/R 
M=10 
MO=M/ 100 
FOR I=1 TO 100 
M1=L2(I)#MO 
M2=13(1)#*MO 
MJ=MO#(BHL2(1)#L2(1)-B-A*Li(I)#L3¢I)) 
Li(l+1)=Licr) +m 
L2(I+1)=L2¢(1)+me 
L3(I+1)=L3¢1)+m3 
LO(I+1)=LO(I) +MO 
NEXT I 
Z1=Z1+1 
IF Zi< 35 THEN 335 
PRINT “AFTER 35 ITERATIONS IN LI@G UNABLE TO MEET REQ’S* 
GOTO 415 
C5=L2(101)-1 
IF H=0 THEN 370 
IF ABS(L3(101)) > .0001 THEN 355 
IF ABS(CS) < .00001 THEN 410 
IF CS > O THEN 365 
H3=H: G6=G3(1): Z=O0: GOTO 140 
HB=H: G7=G3(1): Z=O0: GOTO 140 
IF ABS(L3(101)) > .0001 THEN 380 
IF ABS(CS5) < .COO01 THEN 410 
IF CS > O THEN 350 
G8=G3(1): GOTO 335 
G3=G3(1): GOTO 395 
G3(1)=(G8+G3)/2 
L3(1)=G3(1)2 LOC1)=O: Li¢1)=Git1): Latit)=Lad1) 
GOTO 260 
PRINT "LAMINAR FLOW OVER WEDGE WITH BLOWING" 
PRINT : PRINT : PRINT 
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420 PRINT "B 
425 PRINT “R 
430 PRINT "FG(O) 
435 PRINT "FGDPRIME(O)= ",G3(1),G64,65 

440 PRINT "NI " GO(101),HE,H7 

445 PRINT : PRINT : PRINT 

450 PRINT "NG", "FG", "FGPRIME", "FGDPRIME": PRINT 
455 FOR I=1 TO 101 STEP 10 

460 PRINT GO(I),G1(1),G@¢(1),G3¢1) 

465 NEXT I 

470 PRINT : PRINT : PRINT 

475 PRINT ®*NU". "FL", “FLPRIME", “FLDOPRIME": PRINT 
480 FOR I=1 TO 101 STEP 10 

485 PRINT LO(I),L1(1),LAa¢I), LOI) 

430 NEXT I 

435 PRINT : PRINT : PRINT : PRINT " 7,7 *."END* 


? 8 
Aare 
".G1¢1), "INITIAL LIMITS :" 


1 | | | 
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Figure 6.5 Dimensionless velocity profiles 
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Table 6.1 Comparison Between Wall Shear Stress 
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